\(\vec{PQ}\)
q <- c(5,2)
p <- c(2,3)
pq <- q - p
pq
[1] 3 -1
\(\vec{RS}\)
r <- c(2,-5)
s <- c(1,-2)
rs <- s - r
rs
[1] -1 3
\(proy_{\vec{PQ}}\vec{RS}\)
r <- c(2,-5)
s <- c(1,-2)
rs <- s - r
q <- c(5,2)
p <- c(2,3)
pq <- q - p
proy <- project(rs, pq)
fractions(proy)
[1] -9/5 3/5
Calculamos:
\(u \cdot v\)
u <- c(-2,1,6)
v <- c(2,4,5)
sum(u*v)
[1] 30
\({\|v\|^2}\)
v <- c(2,4,5)
norm(v, type="2")^2
[1] 45
\(w = u - \frac{u \cdot v}{\|v\|^2} v\)
v <- c(2,4,5)
v2 <- norm(v, type="2")^2
u <- c(-2,1,6)
v <- c(2,4,5)
prod_punto_u_v = sum(u*v)
w = u - ((prod_punto_u_v / v2) * v)
w = fractions(w)
sum(w*v)
[1] 0
subspace(as.matrix(w),as.matrix(v))
[1] 1.570796
180*subspace(as.matrix(w),as.matrix(v))/pi
[1] 90
\(W _{\bot } V ?\)
w es ortogonal con v, ya que la múltiplicación entre ellos es igual a cero. Además de que el ángulo que los separa es igual a \(\pi\ /2\)
Con estos puntos: a. Determine si el triángulo \(ABC\) es rectángulo, obtusángulo o acutángulo. b. Determine el perímetro del triángulo \(ABC\) c. Determine el área del triángulo ABC
a <- c(3,0,0)
b <- c(1,0,2)
c <- c(2,3,0)
ab <- b-a
bc <- c-b
ca <- c-a
Respuesta: El triángulo es acutángulo, ya que todos sus ángulos se encuentran entre 0 y 90 grados.
a <- c(3,0,0)
b <- c(1,0,2)
c <- c(2,3,0)
ab <- b-a
bc <- c-b
ca <- a-c
norma_ab <- norm(ab, type="2")
norma_bc <- norm(bc, type="2")
norma_ca <- norm(ca, type="2")
norma_ab + norma_bc + norma_ca
[1] 9.732362
a <- c(3,0,0)
b <- c(1,0,2)
c <- c(2,3,0)
ab <- b-a
bc <- c-b
ca <- a-c
norma_ab <- norm(ab, type="2")
norma_bc <- norm(bc, type="2")
norma_ca <- norm(ca, type="2")
perimetro <- norma_ab + norma_bc + norma_ca
semiperimetro <- perimetro / 2
# Por formula de Héron
area = sqrt(semiperimetro * (semiperimetro - norma_ab) * (semiperimetro - norma_bc) * (semiperimetro - norma_ca))
area
[1] 4.358899
p <- matrix(c(1/2, 1/2, 1/2, 1/2,
(1/sqrt(2)), -(1/sqrt(2)), 0, 0,
(1/sqrt(6)), (1/sqrt(6)), -(2/sqrt(6)), 0,
(1/(2*sqrt(3))), (1/(2*sqrt(3))), (1/(2*sqrt(3))), -(3/(2*sqrt(3)))),
nrow=4, ncol=4, byrow=TRUE)
matriz_p <- fractions(p)
matriz_p
[,1] [,2] [,3] [,4]
[1,] 1/2 1/2 1/2 1/2
[2,] 2378/3363 -5741/8119 0 0
[3,] 19402/47525 19402/47525 -86329/105731 0
[4,] 75658/262087 75658/262087 75658/262087 -489061/564719
p_inversa <- solve(matriz_p)
fractions(p_inversa)
[,1] [,2] [,3] [,4]
[1,] 1/2 2378/3363 19402/47525 75658/262087
[2,] 1/2 -5741/8119 19402/47525 75658/262087
[3,] 1/2 0 -86329/105731 75658/262087
[4,] 1/2 0 0 -489061/564719
p_transpuesta <- t(p)
fractions(p_transpuesta)
[,1] [,2] [,3] [,4]
[1,] 1/2 2378/3363 19402/47525 75658/262087
[2,] 1/2 -5741/8119 19402/47525 75658/262087
[3,] 1/2 0 -86329/105731 75658/262087
[4,] 1/2 0 0 -489061/564719
La matriz P es ortogonal puesto que su inversa y su transpuesta son iguales.
a <- matrix(c(2, -2, -4,
-1, 3, 4,
1, -2, -3),
3, 3, byrow=TRUE)
a
[,1] [,2] [,3]
[1,] 2 -2 -4
[2,] -1 3 4
[3,] 1 -2 -3
a%*%a
[,1] [,2] [,3]
[1,] 2 -2 -4
[2,] -1 3 4
[3,] 1 -2 -3
La matriz A es idempotente puesto que es igual a ella misma al cuadrado.
m <- matrix(c(3/2, -5/2,
2/3, -1/3),
2, 2, byrow=TRUE)
fractions(m)
[,1] [,2]
[1,] 3/2 -5/2
[2,] 2/3 -1/3
m3 <- m%*%m%*%m
m2 <- m%*%m
fx = 6*m3 + 3*m2 - m
fractions(fx)
[,1] [,2]
[1,] -37/6 -55/6
[2,] 22/9 -116/9
a <- matrix(c(1, 2, 3,
2, 5, 7,
-2, -4, -5),
3, 3, byrow=TRUE)
solve(a)
[,1] [,2] [,3]
[1,] 3 -2 -1
[2,] -4 1 -1
[3,] 2 0 1
det(a)
[1] 1
b <- matrix(c(3, -2, -1,
-4, 1, -1,
2, 0, 1),
3, 3, byrow=TRUE)
solve(b)
[,1] [,2] [,3]
[1,] 1 2 3
[2,] 2 5 7
[3,] -2 -4 -5
det(b)
[1] 1
c <- matrix(c(0, 2, 1,
1, 3, -1,
-1, 1, 2),
3, 3, byrow=TRUE)
solve(c)
[,1] [,2] [,3]
[1,] 3.5 -1.5 -2.5
[2,] -0.5 0.5 0.5
[3,] 2.0 -1.0 -1.0
det(c)
[1] 2
d <- matrix(c(3, 6, 9,
2, 5, 1,
1, 1, 8),
3, 3, byrow=TRUE)
det(d)
[1] 0
La última matriz no tiene inversa puesto que el determinante es cero, es decir la matriz es singular o invertible.
Que relación existe entra las matrices que poseen inversas y el valor de su determinante? Sug: revisar la teoría vista en clase.
Autor Brian Duran